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Mahan [J. Math. Phys. 36, 6758 (1995)] has calculated the transmission coefficient and angular 
distribution of particles which enter a thick slab at normal incidence and which diffuse in the slab with 
linear anisotropic, non-absorbing, scattering. Using orthogonality relations derived by McCormick 
& Kuscer [J. Math. Phys. 6, 1939 (1965); 7, 2036 (1966)] for the eigenfunctions of the problem, this 
calculation is generalised to a boundary condition with particle input at arbitrary angles. It is also 
shown how to use the orthogonality relations to relax in a simple way the restriction to a thick slab. 



We consider the equation of radiative transfer with anisotropic scattering in a uniform slab, which occupies the 
space < z < D, together with a boundary fpndition which allows particles to enter the slab through the surface 

z 



at an angle 9 = arccos //q to the normahEf 



l 



f(0,n) = 26{n-Ho) 
f(D,-n) = 



for pi > . (2) 



For thick slabs (D 1), MahanB has presented a solution to this problem which is valid only for (mq = 1. Generalisation 
to arbitrary [j,q is of interest when, for example, the particles which enter the slab come from a point source at finite 
distance, or diffuse before entering the slab. These problems require an integration over the range of incident angles. 
Even for collimated beams, the experimental situation is generally one in which the particles are not normally incident. 
Mahan's method is not readily generalised to solve this problem: his Eq. (78) does not hold when /j,q ^ 1, since then 

The general solution to Eq. ([j]) isU 

l 

f(z,n) = a. + 3j[n-z(l-g 1 )] + J du |^M e -*/« + j(„ _ li )A{v)M L {v)e-'' v 



-l 



where the constants a s and j, and the functions Ml(v) and Mr(v) are to be determined from the boundary conditions. 
The explicit form of the function A(fi) readso 

- -2 = -i 1 - Marctanh^) = - A(/i) , (4) 

where A(/i) is defined by McCormick & Kusceri. To apply the orthogonality relations, it is necessary to rewrite the 
solution in terms of the eigenfunctions used by McCormick & KuscerB: 

^(M) = |P^ + A(i/)*(i/-M), (5) 

which have the property 

<j>- v (fi) = 4> v (-m)< (6) 

Equation ([|) can then be written 
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l l 

f(z, M ) = a s + 3j [fx - z(l - 9l )] + J &vM L (v)Mp)e~ z/v + J AvM^-v)^^)^^^ , (7) 

o o 

where we have absorbed the factor 2/lm into the definition of the functions M r (li) and Ml(li) according to 

M R (p) := - M R (p) , M l (jjl) := - M L ({i) . (8) 

A* A* 

The boundary conditions [Eq. (Q)] then become 

l l 

2<5( M - /x ) = a, + 3j> + J dvM L (u)<l> u (p) + J duM R {-v)ct)^)e- D / v , (9) 

o o 

l l 

= a s - 3j> - 3jD(l - gi) + J &vM L {v)<j)^{p)e- D l» + J &vM R {-v)<j) v (Li) . (10) 

o o 

Defining 



B±{y) ■= \ \m l {v) ± M R {-v) 



and adding and subtracting Eqs. (||) and (|T(]) leads to: 

l l 

S(ji - Mo) = | 3 " S M } =F \jD{l - g x ) + J B±(v)(j> v (fi) du± J B ± {v)e- D ' v ^ v {li) Av . 



(11) 



(12) 



In order apply the orthogonality relations- these equations must be multiplied by a weight function. This function, 
denoted hesse and in McCormick & KuscenO by j(fx), is related, but not identical, to the 7(/i) defined by Mahancl, and 
is given bytl 

^=\xh)-< 0£ " £1 - (13) 

The function X{— fx) can be written in terms of the Ambartsumian function! if)(fi) or the Chandrasekhar _H-function0. 
In the limit c — ► 1 these relationships areElHl3 

Tables of X(-^u), for < /i < 1 are given by Case & ZweifelS; numerical evaluation is straightforward using the 
representationE3 

l 

-c f f cx 2 \ ln(x + fx) 



where c is the albedo for single scattering, equal to unity in the case discussed here. We now multiply Eq. ( |l2| ) by 
7(/i) and integrate over /i from to 1. The integrals over \i can be solved using relations provided by McCormick & 
Kusceru (the numbers above the equals signs in the following refer to the relevant equation numbers): 

7(m) d M = 7o = 1 , (16) 



J 7(mWa* = 7i = ^7o = v = ^o| b=0 = 0.7104, (17) 
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<M M ) 7 (m) d M 6 = 9 , (18) 



o 

J 2 / A , 70 3 v z v 



4 7(1/) 2 



If we denote the extrapolation distance for the Milne problem in the case of isotropic scattering Zo| b _ = 0.7104 by 
simply zq, then, using the above relations, Eq. fll2|) becomes 



Mo 



2 X(- M o) 



{ 3 -} T ^(l- 9l)± /^(^^( 



-v)dv. (20) 



The functions B±(/i) can be calculated by multiplying Eq. ( |l^ ) by </>„< (/x)7(/x) and integrating over /i from to 1. 
Using the orthogonality. relationsEj, one finds inhomogeneous Fredholm equations for B±(fi) which can be solved 
by Neumann iteratiorcJ. In the thick slab approximation, where terms of order e~ D are ignored, these Fredholm 
equations are trivially solved. Equation (pot) for a s and j is then also trivial and independent of B±(n): 



^ f ^ 1 T^(l-.9i). (21) 



2 X(-/xo) I 3jz J T 2 

Once the functions B±(p), and hence Ml(h) and Mn(/j,) have been found, Eq. ( pl| ) provides a s and j and, therefore, 
the density f(fi, z). It is in principle possible to follow this procedure taking into account higher order terms oc e~ D . 
However, the equations become complicated in this case. 

Equations (|l]) enable the transmission coefficient T to be evaluated directly. In terms of the X function we find: 

T = j = — ^ - . (22) 

This result generalises to arbitrary /Uq (0 < Mo < 1) the result of Mahani [Eq. (110)], with which it agrees for fio =3-. 
In the case of isotropic scattering, <?i = 0, Eq. ( p2[ ) is in agreement with the result of McCormick & MendelsonE£l 
[Eq. (35)]. 

Finally, it should be noted that McCormick & Kuscertj have also found orthogonality relations which can be used 
to solve half-space transport problems with higher order anisotropy. 
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